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Abstract.  A new  mechanism  of  vertical  conduction  in  superlattices  with  intentional  disorder  is 
discussed.  We  show  that  at  low  temperatures  the  conductance  of  these  structures  can  be  mostly 
determined  by  phonon-assisted  transitions  between  the  second-nearest  wells  in  the  vicinity  of  the 
well  with  the  highest  size-quantization  level.  Such  transitions  involve  virtual  transitions  to  the  state 
of  the  intermediate  well  and  are  characterized  by  a low  activation  energy  and  specific  dependence 
of  the  vertical  conductance  on  the  scale  of  disorder. 


1 Introduction 

Superlattices  with  intentional  disorder  (SLID)  and  multiple  quantum  well  structures,  in 
which  disorder  was  introduced  by  random  controlled  variations  of  the  well  widths  in  the 
process  of  structure  deposition,  were  first  discussed  in  | ] and  were  experimentally  realized 
in  I . ].  Vertical  conduction  (in  the  direction  of  the  SLID  growth  axis)  in  such  structures 
was  studied  by  optical  methods,  in  particular,  stationary  and  picosecond  luminescence  spec- 
troscopy, and  also  by  direct  measurements  of  the  vertical  conductance.  The  conductance 
measurements  for  Si-doped  GaAs/GaAlAs  SLID  revealed  some  unusual  features  [ , ]. 
Thus  at  low  temperatures  the  temperature  dependence  of  the  conductance  was  practically 
nonactivated  (quasimetallic)  even  for  structures  with  large  disorder,  namely,  when  the  width 
of  the  distribution  of  size  quantization  levels  exceeded  the  estimated  miniband  width.  Even 
though  the  Coulomb  fields  arising  from  the  electron  redistribution  between  the  wells  can 
give  rise  to  an  appreciable  suppression  of  disorder  (narrowing  of  the  level  distribution)  [ ], 
the  explanation  of  the  relatively  weak  temperature  dependence  of  the  conductance  still 
remains  problematic. 

The  wave  functions  of  electronic  states  in  SLID  can  be  written  in  the  form  t/^k(|  — 
NU\(z)  exp  (7k|jp),  where  N is  the  normalization  factor,  z is  the  coordinate  in  the  SLID 
growth  direction,  p is  the  in-plane  position  vector,  Ux  (z)  is  the  eigenstate  corresponding  to 
the  solution  of  the  one-dimensional  problem  with  the  potential  V (z)  = J2n  Vi  (z)  describing 
the  modulation  of  the  conduction  band  edge,  and  Vn  (z)  is  the  potential  of  the  n th  well  (we 
set  Vn (z)  = 0 in  the  barrier  regions).  The  functions  Ux(z)  are  localized  and  for  small 
overlap  of  wave  functions  of  the  neighboring  wells,  we  can  use  the  basis  of  “atomic-like” 
wave  functions  localized  at  the  corresponding  wells.  Since  for  the  structures  of  the  type 
studied  in  [ ] the  contributions  >from  higher  subbands  are  negligible,  we  can  take  only  the 
lowest  subband  into  account  and  write  a —n.  For  SLID  the  vertical  conduction  is  usually 
controlled  by  phonon-assisted  tunneling  between  neighboring  wells.  As  in  the  standard 
hopping  theory  (e.g.,  see  [ , ]),  the  low-field  problem  can  be  reduced  to  the  equivalent 
resistance  network,  with  resistances  expressed  in  terms  of  resulting  transition  rates  between 
the  wells  (the  difference  from  the  standard  hopping  problem  is  that  we  have  to  sum  over 
the  initial  and  final  states  of  the  wells).  Thus,  for  the  resistance  Rnni  connecting  the  wells 
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n and  n'  we  have 

R~n'  = (e2/ kT ) E ^nkn.n'k'/nkjCl  ~ A'k'b  (D 

k||,k'| 

where  ky  is  the  two-dimensional  momentum  in  the  well  plane,  W„k|| is  the  probability 

of  phonon-assisted  transitions  from  the  state  kj  to  the  state  «k  and  /„  kB  is  the  equilibrium 
average  occupation  number  of  the  state  uky.  If  the  higher  level  (say,  E„)  lies  above  the 
Fermi  level  n,  then  (for  En  — (t»  kT)  we  have  Rnnr  — Ro  exp  {—{En  — n)/kT},  where 
the  preexponential  factor  Ro  depends  on  the  overlap  of  the  wave  functions  of  neighboring 
wells  and  only  weakly  depends  on  energy  and  temperature.  In  the  nearest-neighbor  ap- 
proximation, the  network  is  quasi-one-dimensional,  and  its  resistance  is  given  by  the  sum 
of  series  resistances  connecting  of  the  neighboring  wells  of  the  chain.  Clearly,  if  the  scatter 
of  the  energy  levels  is  greater  than  the  miniband  width  and  kT,  then  the  total  resistance  is 
determined  by  a critical  well  (or  by  a small  number  of  critical  wells)  with  adjoining  largest 
resistances.  These  resistances  exponentially  depend  on  temperature  and  it  follows  that  the 
temperature  dependence  of  the  total  vertical  conductance  is  activated,  with  the  activation 
energy  determined  by  the  position  of  the  highest  energy  levels  relative  to  the  Fermi  level. 
However,  with  lowering  temperature  the  vertical  conductance  can  be  controlled  by  parallel 
shunt  resistances  corresponding  to  second-neighbor  phonon-assisted  transitions.  In  fact, 
the  critical  resistances  exponentially  increase  and  the  level  separation  for  transitions  be- 
tween wells  lying  on  different  sides  of  the  critical  one  (and  hence  the  activation  energy) 
is  typically  smaller  for  these  transitions.  We  shall  see  that  the  most  probable  transitions 
involve  virtual  intermediate- well  states. 

2 Probability  of  second-neighbor  phonon-assisted  transitions  via  virtual  states 

Let  us  consider  a three-well  configuration  representing  the  critical  region  consisting  of  the 
critical  well  2 with  a high  level  En  lying  above  the  Fermi  level  and  two  adjacent  wells  1 and  3 
with  lower  levels.  Next,  let  un  (z)  and  En  be  the  wave  functions  and  energies  corresponding 
to  the  solution  of  the  one-dimensional  problem  with  a single  nth  well  (we  assume  that 
En  > E]  > E i).  Taking  the  overlap  into  account,  we  can  construct  the  hybridized 
wave  functions  corresponding  to  the  solution  of  the  three-well  problem.  Thus,  for  level 
separations  exceeding  the  transfer  integrals  tmn  — f dz  um(z)Vn(z)  un(z),  where  Vn(z)  — 
J2n'^n  Vb'(z),  the  lowest  energy  state  is  mostly  localized  at  well  1 and  the  corresponding 
wave  function  is 

U\{z)  - Ni{ui(z)  + cnU2(z)},  (2) 

where  N i is  the  normalization  factor  and  c 12  = t\i/t E2  — E\).  Here  we  have  omitted 
the  transfer  integrals  between  wells  1 and  3 and  as  usually  neglected  the  nonorthogonality 
integrals  with  the  functions  un  (z).  Similarly,  we  can  write  out  the  expressions  for  the  wave 
functions  Un(z)  and  Un,(z)  localized  at  the  wells  2 and  3. 

The  phonon-assisted  transition  probability,  which  involves  the  spatial  displacement 
between  wells  1 and  3 is  expressed  in  terms  of  the  matrix  element 

'3k| . ik;  = f dzUiMH^U^z) 

— NiN3{cn  j dzu2(z)He'Jun,(z)  + C32  J dzu\ (z)He'lJ  uiiz)},  (3) 
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where  He  — J dp  exp  { ? (k  n — k^)p}Hepi,  and  He  pi,  is  the  Hamiltonian  of  the  electron- 
phonon  interaction. 

Expression  (3)  describes  the  amplitude  of  the  transition  between  wells  1 and  3 as  the  sum 
of  contributions  of  the  two  channels,  the  first  one  corresponding  to  the  succession  of  virtual 
tunneling  from  1 to  2 succeeded  by  phonon-assisted  transition  from  2 to  3 and  another  one 
corresponding  to  phonon-assisted  transition  from  1 to  2 succeeded  by  tunneling  from  2 to  3. 
Note  that  for  our  system  the  probability  of  direct  phonon-assisted  transition  from  1 to  3 is 
small  compared  to  that  of  the  process  described  by  Eq.  (3).  In  fact,  the  integrals  involving  the 
functions  u\ (z)  and  113(7.)  are  proportional  to  the  overlap  factor  exp  {— or(2u>  + £)},  where 
a is  the  inverse  wave  function  decay  length  in  the  barrier  region,  w is  the  barrier  width 
and  L is  the  well  width.  On  the  other  hand,  the  product  of  the  overlap  factors  appearing 
in  Eq.  (3),  exp  (—2 aw),  is  much  greater  so  that  for  exp  (—o'/.)  <?c  1 the  probability  of  the 
direct  transition  from  1 to  3 is  negligible.  This  is  in  contrast  to  the  conventional  situation 
in  the  problem  of  hopping  between  localized  states,  where  usually  the  probability  of  direct 
phonon-assisted  transition  between  distant  sites  is  much  larger  than  that  for  transitions 
involving  virtual  states. 

The  shunt  resistance  R\3  for  the  processes  involving  virtual  states  is  defined  by  ( I),  where 
the  transition  probabilities  are  expessed  in  terms  of  the  modulus  of  the  matrix  element  (3) 
squared.  Provided  that  (£3  — £|)  > kT,  we  obtain  that  the  conductivity  is  activated  with 
the  activation  energy  £3  — \i,  which  is  smaller  than  the  activation  energy  Ei  — n for  the 
conduction  mode  with  nearest-well  transitions  in  the  critical  region.  Another  difference 
between  the  nearest-well  and  distant-well  conduction  modes  is  a different  scaling  behavior 
as  the  disorder  energy  (the  width  of  the  level  distribution)  is  varied.  As  follows  >from 
the  corresponding  dependence  of  | e 1 2 1 2 , the  preexponential  factor  of  the  conductance  is 
inversely  proportional  to  the  disorder  energy  squared  whereas  the  activation  energy  scales 
as  the  disorder  energy.  Note  that  the  usual  experimentally  studied  structures  with  the 
number  of  the  wells  of  the  order  100  [ ] are  mesoscopic  and  can  have  widely  fluctuating 
conductances  for  various  realizations  of  the  size  quantized  level  distribution.  The  scaling 
behavior  of  the  conductance  can,  however,  be  checked  for  these  structures  if  one  uses  series 
of  structures  similar  to  those  described  in  [ ] with  different  disorder  energies  for  the  same 
random  realization. 

3 Discussion 

The  scaling  behavior  of  the  preexponential  factor  of  the  conductance  is  expected  to  become 
evident  in  the  case  of  small  activation  energies.  The  activation  energy  can  be  small  if  the 
energies  £1,  £3  for  the  states  localized  on  different  sides  of  the  critical  well  lie  below  the 
Fermi  level.  Actually,  the  levels  £j , £3  can  correspond  to  hybridized  states  of  several  wells 
(clusters)  rather  than  to  individual  wells  adjacent  to  the  critical  one.  If  the  miniband  width 
for  the  corresponding  regular  SL  (or  the  transfer  energy)  is  not  too  small  compared  with  the 
disorder  energy,  then  the  average  number  of  sites  in  the  clusters  is  sufficiently  large;  then 
for  the  most  probable  realization  of  disorder  the  lowest  energies  of  hybridized  states  of  the 
clusters  adjoining  the  critical  well  lie  below  the  Fermi  level.  In  this  case  the  nonactivated 
temperature  dependence  of  the  vertical  conductance  is  expected. 
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